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O ■ ABSTRACT 
T— I \ 

\ We examine the necessity of resolving the Jeans mass in hydrodynamical simula- 

^r^ ' tions of the hierarchical formation of cosmological structures. We consider a standard 

' two-component fluid and use smoothed particle hydrodynamics to model the baryonic 

,—1 ■ component. It is found that resolution of the Jeans mass is not necessary to extract 

^ ■ bulk properties of bound objects that are independent of resolution, provided that the 

. objects have undergone a number of merger events. The degree of merging must be 

^ ! sufficient to populate the objects with substructure of similar total mass. After this 
criterion is met, neither the density profiles of structures nor their density-temperature 

O ' distributions depend appreciably on the resolution of the simulation. Resolution of the 

' Jeans mass by the imposition of a minimum temperature is not found to adversely af- 

■ feet the state of the matter in the final clusters. The baryon concentration profiles still 

■ indicate a deficit of gas with respect to the dark matter in the interiors of structures. 

^: 

. Subject headings: galaxies: clusters: general — hydrodynamics — large-scale structure 

of universe — methods: numerical 



1. INTRODUCTION 

Numerical simulations have proved useful for understanding the history and dynamics of cos- 
mological structure (Davis, Efstathiou, Prenk, & White 1992; Katz 1992; Cen, Miralda-Escude, 
Ostriker, & Ranch 1994; Evrard, Metzler, & Navarro 1996). Our ability to model collisionless 



Present address; University of Maryland, Baltimore County, Department of Physics , Baltimore, MD, 21250 

^Present address: McMaster University, Department of Physics and Astronomy, Hamilton, Ontario, L8S 4M1, 
Canada 



- 2 - 



systems has matured to the extent that bulk properties of the dark matter can be predicted with 
increasing confidence, although important exceptions exist, such as the merger rates of substructure 
as well as the inner profiles of dark matter halos (Moore, Governato, Quinn, Stadel, & Lake 1998, 
see, for example,). 

Combining collisionless N-body techniques with a hydrodynamical method to follow the bary- 
onic component opens the way to a detailed understanding of cosmic structure formation in a 
variety of models. Modeling the baryonic component removes the need to make assumptions about 
the relationship between the baryonic material and the underlying collisionless matter. Such tech- 
niques have proved to be both popular and powerful (Hernquist &; Katz 1989; Evrard 1990; Katz & 
Gunn 1991; Cen 1992; Navarro & White 1993; Steinmetz & Miiller 1993; Ryu, Ostriker, Kang, & 
Cen 1993; Bryan, Cen, Norman, Ostriker, &; Stone 1994; Gncdin 1995; Steinmetz 1996). A recent 
comparison of several hydrodynamic codes (Frenk et al. 1998) has shown that we may also begin 
to have confidence in the overall properties of the baryonic distribution predicted by these codes. 
However, several differences which are apparent both between codes and between results of varying 
resolution suggest that convergence has yet to be established conclusively. We address here the 
particular issue of convergence with resolution in hydrodynamic simulations in which objects are 
built up through hierarchical merging. We do not consider radiative cooling in this work. 

In the standard model of cosmological structure formation, larger structures are formed via 
the amalgamation of smaller structures formed at earlier times in a hierarchical fashion. In any 
simulation which models the formation of such structures in this hierarchical manner, there will 
be an abundance of objects formed at the resolution limit of the simulation. Numerical effects will 
be a particular problem for these objects. The propagation of these effects into the larger objects 
subsequently formed requires characterization both to determine the reliability of the physical 
model as well as to ensure convergence among simulations performed at various resolutions. 

Recently Owen & Villumsen (1997) has argued that it is necessary to resolve the Jeans mass 
in a simulation if convergence is to be attained with increasing resolution. This requires that a 
minimum temperature be imposed in the simulation such that the gas in the first bound objects — 
which form at the resolution limit in a standard simulation of hierarchical clustering — is at a 
sufficient temperature to support itself via pressure in the gravitational potential well of the first 
dark matter structures. They argue that if the gas temperature is initially much less than the virial 
temperature of the first halos (this is a boundary condition frequently adopted in cosmological 
simulations), the cold gas falling into the first potential wells will suffer spuriously large amounts 
of shock heating of the gas. This "shocking" will be greater for simulations with lower resolution 
since the scale over which the shocks occur will be larger. Conversely, if the initial generation of 
objects is pressure-supported, the interaction of the baryons in merging halos as the hierarchical 
structure formation proceeds will lead only to weak shocks. (This result follows since the velocity 
with which lumps merge is similar to the internal velocity dispersion of the individual lumps in a 
typical hierarchy.) Owen and Villumsen suggests that the difference in the amount of shock heating 
at early stages will persist through to the final object and prevent convergence of numerical results 
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unless a minimum temperature is imposed. 

Owen k. Villumsen finds, in scale-independent simulations, that without the imposition of a 
minimum temperature to pressure support resolution-limited structures, the gas is found to move 
en masse toward a state of higher density and lower temperature as the resolution of the simulation 
increases. This behavior is not to be expected during the formation of hydrostatic structures. In 
the standard model, we expect the virial temperature of a structure to be broadly determined by 
its mass, following the relationship T oc M^/^. Since the bulk of the gas resides in the halos of these 
structures, the temperature of the bulk of the gas should consequently be determined by the dark 
matter distribution, a distribution whose convergence among differing resolutions is confirmed in 
the simulations of Owen & Villumsen. 

There are other issues of concern raised by the results of the simulations described in Owen 
&; Villumsen (1997). First, the results suggest paradoxically that the injection of thermal energy 
into a system at an early time can lead to a cooler system at a later time. The authors claim 
that this is achieved by reducing the amount of shock-heating as the first, unresolved, objects 
collapse. However, the injection of thermal energy into the system via the imposition of a minimum 
temperature pre-heats the gas in a manner which is essentially equivalent to shock-heating the gas 
in the first objects. 

Another concern is the concentration of baryons relative to the dark matter in the cores 
of structures reported in Owen &; Villumsen. In three-dimensional simulations without cooling, 
spanning a wide range of resolutions, others have found extended distributions of baryonic material, 
relative to the dark matter (Evrard 1990; Thomas & Couchman 1992; Cen & Ostriker 1993; Kang, 
Cen, Ostriker, k. Ryu 1994; Metzler & Evrard 1994; Pearce, Thomas, & Couchman 1994; Navarro, 
Prenk, & White 1995; Anninos & Norman 1996; Lubin, Cen, BahcaU, & Ostriker 1996; Pildis, 
Evrard, & Bergman 1996). Pearce et al. (1994) explains the phenomenon as a result of the merging 
process in which gas is shocked, permanently removing energy from the dark matter component 
and passing it to the gas. 

A possible explanation for the conflicting results regarding the baryon concentration may be 
that the simulations of Owen & Villumsen were done using a 2-D code. This allowed them to 
achieve high spatial resolutions but may have introduced phenomena unique to two-dimensional 
geometries. 

1.1. Jeans Mass in a Two-component Fluid 

In a self-gravitating baryonic fluid, the Jeans mass is easily obtained, by a perturbatitive anal- 
ysis, as a sharp threshold between acoustic oscillation and collapse, depending on the density and 
temperature of the fluid. In a two-component fluid comprised of baryons and coUisionless cold dark 
matter, the situation is more complicated. The low initial velocity dispersion of cold dark matter 
leads to instability, or growth, in this component on all scales. The relevant question in this case for 
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the baryonic fluid is on what scales do the two components remain coupled adiabatically? On small 
scales, gas pressure will be sufficient to oppose the gas falling into the small dark matter potential 
wells, whereas on large scales the initial gas temperature will be lower than the effective virial 
temperature of the forming dark matter halo and the gas will increase in overdensity with the dark 
matter. A perturbative analysis of the coupled two-component fluid is relatively straightforward 
but not very illuminating for our purposes. In particular it can be shown that there is not a single 
Jeans mass in the sense described above (see Peacock 1999, for example). It is sufficient for our 
purposes to define the Jeans mass as that which is just sufficient to adiabatically compress the gas 
to an overdensity which is the same as that of the dark matter. 

With this assumption we can approximate the Jeans mass of the gas for a fluid composed of 
collisionless dark matter of density, poM, and collisional baryonic material (gas) of density, pg, and 
temperature, T, by 

2 \G pM P^ruH J 
where the total mass density is given by pM = Pdm + Pg- 



1.2. The Minimum Temperature for Pressure Support 

The effective mass resolution of a simulation is usually fixed. For simulations such as those 
described here which use smoothed particle hydrodynamics (SPH) to model the gas forces, this 
can be taken to be the total mass of the number of particles, Nsph, over which local quantities 
are averaged, presuming all gas particles have equal mass. The number of particles is typically on 
the order of 30. By associating the mass of this number of particles to the Jeans mass (Eq. 1) and 
assuming a total mass density, pM in Eq. 1, equal to the mean density of the simulation volume, 
Pjy^, the temperature required to support gas in dark matter halos of this size is 

_ / 2NsPHmg \ ^'^ fuuH 
\ Pg J k 

where m,g is the mass per gas particle. If the gas has a minimum temperature set by this value, then 
we can say that the simulation is resolving the Jeans mass. This assumption is flawed in a number 
of ways. Since local density is not in general equal to the mean in the simulation, the minimum 
temperature should also be a function of position. However, the minimum temperature should 
vary with the local density as T^i„ — Pmi which is not a particularly sensitive function, especially 
in those regions in which the first objects are forming. More significant is the decoupling of the 
components of the fluid, discussed previously. This is likely to occur for the first objects formed 
since the dark matter halos collapse regardless of whether or not the gas is pressure-supported. 
For those situations in which, locally, poM ^ Pg-, enhancement of the dark matter density without 
a similar increase in the gas density requires the temperature to rise by an equivalent fraction 
in order to maintain a constant Jeans mass for the gas. However, during the formation of these 
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dark matter structures, the gas will be allowed to gradually respond, adiabatically, to the increasing 

gravitational potential, partially or fully offsetting the need for an increased minimum temperature. 
This would not occur without the initial minimum temperature since the gas in that case would 
initially collapse along with the dark matter. 

1.3. The Hypothesis 

In this paper, we propose that the effect on the state of the gas of poorly resolving the first 
baryonic structures does not propagate into larger structures formed through hierarchical merger. 
The properties of the gas in the larger halos will be determined primarily by the virial properties 
of the structures which are dominated by the distribution of the collisionless component. 

To test this hypothesis, we performed a series of seven simulations of hierarchical structure 
formation. The simulations were done at three different resolutions. For two of these simulations, 
the temperature of the gas content was prevented from being less than the temperature required to 
maintain pressure support for the smallest resolvable object; that is, the Jeans mass was resolved. 

Comparison of the data simulated for this study was made using both the total sample of gas 
in each simulation volume as well as an individual wcll-rcsolvcd structure. Using an individual 
cluster permits examination of those effects that are frequently the focus of interest in structure 
formation simulations; the void particles in simulations of this sort are usually ignored due to a 
lack of resolution in this regime. 

We examine the entropy of the gas, which is a measure of the degree of shocking. Also examined 
is the density distribution of the gas to determine the extent of any bulk movements of the gas to 
higher densities as resolution is increased. We also explore mean bulk properties of the structures: 
in particular, the mean dark matter and gas density profiles and the mean baryon fraction profile. 
The mean gas density and baryon fraction profiles are potentially sensitive to a lack of convergence 
as well as being of fundamental interest. 

The outline of the paper is as follows. The simulations are described in Sec. 2. The results of 
an attempt to recreate the results of Owen &: Villumscn but using a three dimensional geometry 
are given in Sec. 3.1. The impact of resolving the Jeans mass on the final state of the gas in galaxy 
clusters is addressed in Sec. 3.2. The results are discussed in Sec. 4. Hereafter, we refer to Owen 
& Villumsen (1997) as OV97. 

Though the simulations are scale invariant, the results on occasion are given in physical units. 
This is simply a matter of convenience for the purpose of comparison. 
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2. 



THE SIMULATIONS 



The N-body code with SPH, hydra (Couchman, Thomas, & Pearce 1995), was used for all 
simulations. HYDRA calculates the gravitational forces using an adaptive mesh for the large-scale 
complemented with particle-particle calculations for the short range forces. Smoothed particle 
hydrodynamics (SPH) is used for the hydrodynamics. 

Seven simulations were performed. Five used the original code, and two used a variant of 
HYDRA which imposes a minimum temperature on the gas particles. The details of the simulations 
are outlined in Table 1. Given, there, are the number of each type of particle and the gravitational 
softening length, e. The simulations that were evolved to the final expansion factor, a = af, were 
performed twice, once with and once without a minimum temperature, whereas the first trio were 
evolved only once, without a minimum temperature. Initial density perturbations axe the same for 
all resolutions within a group, albeit limited by the Nyquist frequency. 

The minimum temperature is set to support a mass of particles corresponding to the effective 
mass resolution of the simulation as described in the introduction. At each epoch of the simulation, 
structures composed of Nsph gas particles with a local density equal to the mean gas density of 
the volume, p^, and individual masses of nig are supported against collapse by the imposition of a 
minimum temperature, Tmin- The minimum temperature is derived by simply equating the Jeans 
mass, Mj, to the structure mass NspH^g, as is done in Sec. 1.2 to derive Eq. 2. 

During the two simulations in which a minimum temperature was imposed, Tmin was recalcu- 
lated at each epoch. Since the minimum temperature scales with the mean density as Tmin P^^^) 
then it scales with the expansion factor, a, as T^m oc a~^. Normally, adiabatic cooling for a uni- 
form monotonic gas would maintain the relation, T oc a~^. Consequently, the imposition of the 
minimum temperature has the effect of continuously injecting thermal energy into the fraction of 
the gas whose adiabatic expansion would bring its temperature below the minimum temperature. 

The initial conditions for the two-dimensional simulations of OV97 have a density perturbation 
spectrum with a spectral index equivalent to the three-dimensional n = —1 form used in the 
simulations presented here. Because the simulations are scale free, the spatial variance of the 
density fluctuations at the end of the simulation in the larger volume corresponds to an earlier 
epoch in the smaller volume, after scaling by the box size. 

The amplitude of linear density fluctuations, as measured by the power spectrum of the density 
field, P{k), relates to the wavenumber of the fluctuation, k, measured in comoving coordinates, and 
the time, t, via 



where ctl is a normalization factor corresponding to the assumed value of the variance, a, of the 
present density field on the scale of L. 




(3) 



Although the simulations discussed in OV97 and those presented here are scale invariant, it 
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is useful to set a physical size to the simulation volume. Indeed, this may be a benefit to those 
readers familiar with the scale of cosmological structures. Further, this is the only practical way we 
have of determining the degree of clustering in the simulations of OV97. OV97 uses a volume with 
a side of length 128/i~^ Mpc , h = 0.5 (Owen &: Villumsen 1996). Since we use a volume with sides 
40/i~^ Mpc and h = 0.65, this gives a ratio of scales of 4.16 for their simulation scale compared to 
the simulations described here. For the normalization factor, they have erg = 0.6 whereas we used 
(Tg = 0.935. This implies an expansion factor of 0.24a/. That is, the strTictures in the simulations 
presented here, scaled by the box size, will have grown in similar magnitude to those of OV97 after 
an expansion of 0.24 of the final expansion factor of the simulation. 

In this study, the data have been compared at the slightly earlier expansion factors of 0.17a/ 
and 0.22a f to ensure that structures have not undergone more than a couple of merger events. As 
a check, a volume with the same cosmological parameters as described in OV97 but with the same 
perturbation waves as found in our smaller boxes was evolved to an expansion factor, a/. The data 
from the the 40^^^ Mpc volumes at an expansion factor of 0.22a/ are visibly similar to the final 
output of this test volume, supporting the legitimacy of scaling the numerical results. 

It was feasible, computationally, to evolve the volume to the present only at the resolutions of 
2 X 64^ particles and less. Hence, this study will concentrate on a comparison of the simulations of 
resolution corresponding to this number of particles with simulations of a resolution appropriate to 
2 X 32^ particles. These resolutions will be referred to henceforth as the 64^ and 32"^ simulations. 
Each was evolved with and without a minimum temperature. For the earlier epoch corresponding 
to the expansion factor of 0.17a/, the volume could also be evolved at the resolution provided 
by 2 X 128^ particles. Thus a set of three simulations was performed using the same density 
perturbations (limited by the Nyquist frequency) at resolutions corresponding to 128^, 64^, and 
32^. 

Cooling was neglected to maintain the scale invariance. In all simulations, $7 = 1, 0,dm = 0.9, 
and Q,g = 0.1. 

The initial density perturbations were established by displacing the particle positions from a 
uniform cubic grid using the Zel'dovich (1970) approximation for growth of density perturbations 
in the linear regime in the standard way. 

The initial redshift for the 32^ and 64'^ simulations is Zinit = 75. This was chosen to keep the 
maximum displacement incurred during the establishment of the density field to less than 1/2 of the 
initial grid spacing. This ensures that the correct non-linear field will be attained at later times. 
The same criterion was used to set the initial redshift for the 128'^ simulation for which Zinu = 150. 
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3. RESULTS 

3.1. The State of the Baryons at an Early Epoch 

We first attempt to reproduce the results of OV97 by examining the state of the baryons in 
our simulations at a comparable epoch with that of the final epoch of the simulations of OV97. The 
smaller physical scale of the simulations presented here requires that we examine our data before 
the nominal end of the simulation in order that a comparable level of clustering has developed 
to that present at the final epoch of the simulations of OV97. The required expansion factor is 
a = 0.22a f, where a/ is the final expansion factor of the simulations. 

Following the analysis given by OV97, we look at the distribution of the gas particles in the 
p — T plane as well as the density distributions for which the density is resampled down to the 
resolution of the lowest resolution runs as described below. The position within the p — T planes 
gives an indication of the degree of shocking, with the more highly shocked gas (gas with greater 
entropy as indicated by Tp~'^, 7 = 2/3) lying to the upper left (high temperature, low density). 

In order to allow a direct comparison of the densities in simulations of different resolution a den- 
sity resampling technique was employed. The (resolution-degrading) resampling was accomplished 
by smoothing over equal masses for each particle using an SPH-like density estimate. Consequently, 
there are approximately eight times as many particles in each estimate for the 64^ data as in the 
32^ data, but the smoothing radii are approximately the same. This resolution degradation is 
performed for the high resolution simulations (64^) when comparison is made with the lowest reso- 
lution simulation (32^). This techniqiic was unfortunately not possible with the 128^ data-set due 
to the computational effort of summing over the more than 3000 neighbor particles of each gas 
particle. 

It is clear that for density profiles without cores, higher mass and spatial resolutions will lead 
to more gas at a higher density even if the simulations can be said to have converged to the true 
solution within resolution limits. If the densities are averaged within volumes encompassing a mass 
which is resolved in the lowest resolution simulation, then these masses should not vary significantly. 

The simulations of OV97 produce more cool, dense gas as the resolution increases. They 
interpret this to mean that in the higher resolution simulations, the gas in the cores of structures is 
not as strongly shocked. To search for this effect, three sets of data were evolved at the resolutions 
provided by 2 x 32^, 2 x 64^, and 2 x 128^ particles, half in each phase of gas and dark matter, 
as described in Sec. 2. Each set has the same initial perturbations, limited appropriately in the 
high spatial frequency domain by the Nyquist frequency. These data sets were evolved to an early 
epoch for which the amplitude of density fiuctuations has grown to a level comparable to those in 
the data presented by OV97. 

The evolved data were resampled, in the case of the 64'^ and 128^ sets, to the same mass 
resolution as the 32^ set by randomly selecting 32^ particles from each phase. This avoids the 
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problem of performing the density estimate over the very large number of neighbor particles required 

for the 128^ data set, at the expense of an increase in shot noise. Local SPH density estimates for 
the sampled particles in these resampled sets were determined allowing comparisons of the data in 
p-T space. 

The p-T distributions for all the gas in each of these three simulations are illustrated in Fig. 1 
(top panels). This is comparable to the top panels of Fig. 8 of OV97. It is clear that there is no 
trend for the hulk of the mass to move toward a cooler and denser state as the resolution increases, 
in contrast to the results of OV97. Since the bulk of this gas is in the halos of bound structures, 
the argument given in the introduction is supported; the entropy of the gas is set by the virial 
properties of the halo which are dictated by the dark matter. There is a trend for the void particles 
(found to the lower left in the p-T plane) to increase in temperature as the resolution increases. 
There is also a trend for the distributions to become more diffuse, particularly toward the denser 
and cooler state. 

The lower set of panels in Fig. 1 give the distribution for particles in a cluster common to all 
three sets of data. The cluster was selected simply for being one of the larger clusters of particles 
to have formed by this early stage. The cluster is barely resolved (~ 100 particles) in the 32^ 
simulation, but is well resolved in the 128^ simulation. Here it is evident that the diffusion of 
the distributions to a denser and cooler state is due to a small population of particles, and not a 
bulk movement. Upon examination of the spatial distribution of this population, these particles 
are found to be located in substructure that is present only in the simulations with greater mass 
resolution. These clumps are not as shocked as the particles in the halos of the substructure. This 
is consistent with the phenomenon described by OV97 in nature, though not in extent since it is 
not a bulk motion, but simply a new population of particles. 

As the mass resolution increases, the first bound haloes form at earlier epochs at higher densi- 
ties leading to increased central densities in the final dark matter halos. This is illustrated in Fig. 2 
which compares the gas and dark matter mean density profiles for resolved clusters in the simula- 
tions, albeit at a later epoch. The dark matter density is double in the 64'^ simulation compard with 
the 32'^ simulation in the central regions (close to the resolution limit of the 32^ simulation) . This 
increases the degree of shocking of the halo particles more than offsetting the decreased amount 
of shocking in the cores of substructure. The mean temperature of this structure thus increases 
with increasing resolution: in physical units it is 2.9 x 10^ K , 3.2 x 10^ K , 4.5 x 10^ K for the 32^, 
64^, and 128'^ simulations, respectively. The total amount of shocking on the whole also increases, 
but not to such a degree. The amount of shocking, as measured by the mass-weighted mean of 

2/ 

the entropy parameter, < T/pg >m, for an individual cluster at the three resolutions evolved to 
a = 0.17a/ is given in Table 2. The values not resolution-degraded are given in the second column 
which reveal a slight increase (10% for every doubling in resolution) in the entropy as the resolu- 
tion increases. If the densities are recalculated for the resampled data with the intent of comparing 
similar effective resolutions, then the degree of shocking decreases quite markedly with increasing 
resolution (Table 2, third column). This occurs simply because the temperatures of the particles 
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are not changing while the resampled densities of the particles in the densest regions decrease due 

to the use of a larger smoothing radius. It would be erroneous to conclude that the gas is actually 
less shocked, since the temperature of the halo gas particles and the depth of the potential well are 
so intimately related. 

The gas density distributions for both the 32'^ and 64^ data were found at an expansion factor of 
0.22a f. The densities of the 64'^ data-set were resampled as described previously. The distribution 
of gas density at this early epoch (Fig. 3, left hand panels, dashed lines) is only weakly dependent 
on resolution, becoming somewhat less peaked at 10~^^ g cm at the higher resolutions. No bulk 
motion of the gas to higher densities is observed, unlike the simulations described in OV97. 

Since no bulk motion of the gas to higher densities and lower degrees of shocking is observed, 
we might predict that the imposition of a minimum temperature will not significantly affect these 
results. The distribution of the gas within the p — T plane at the early epoch of a = 0.22a f is 
illustrated in Fig. 4. Again, the bulk of the gas is found in the same state in both simulations 
without a minimum temperature (Fig. 4, left panels), with a diffusion toward a state of lower 
entropy. For those simulations in which a minimum temperature is imposed (Fig. 4, right panels), 
the gas is displaced within the p — T plane in those regions below T^i„ (for obvious reasons) and 
below the locus that uniform gas would trace if it were allowed to cool to the minimum temperature 
appropriate at each epoch. In order for gas to be found below this locus it must be initially held at 
a high temperature and density but with little or no increase in entropy and subsequently allowed 
to cool adiabatically. The bulk of the gas is above these boundaries. The gas that is above these 
boundaries is consistently located in similar proportions to the gas found in the same locations in 
the p — T plane in the simulations without a minimum temperature. There is no gas in the lower 
resolution simulations which fits these criteria, however. 

The resampled density distributions for the simulations with a minimum temperature, illus- 
trated in Fig. 3 (right panels, dashed lines), confirm that the minimum temperature does modify 
the distributions but docs not provide convergence. This figure is comparable to Fig. 6d in OV97 
which indicates a bulk motion to higher densities as the resolution is increased, even when the 
data are resampled to the scale of the lowest resolution simulation. In our simulations, for the gas 
denser than 10~^^g cm~^, there is actually a decrease in the agreement between the lower and 
higher resolution data when a minimum temperature is enforced. 

The decrease in the degree of shocking reported by OV97 as the resolution of a hydrodynamic 
simulation is increased is confirmed in the cores of clusters of particles. However, the bulk move- 
ment of the gas particles to a less shocked state is not observed in our simulations. Indeed, the 

mean degree of shocking of cluster particles is essentially conserved, with the mean temperature of 
particles actually increasing with resolution. The displacement in density distributions reported by 
OV97 is also not observed, nor is any sort of convergence provided by the imposition of a minimum 
temperature. 
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3.2. The State of the Baryons at a Late Epoch 

After many mergers, it is possible that the gas in the initial structures described above will 
have 'forgotten' the degree of shocking it had at early times. It may have been sufficiently shocked 
by subsequent mergers that early heating, via strong shocks or otherwise, would no longer be 
significant to the total thermal content and distribution. To examine this hypothesis, the 32^ and 
64^ sets of simulations used in Sec. 3.1 were evolved to a further state which involved many levels 
of merging. It is not presently feasible to evolve a 128^ to such a state with the hardware on hand 
in a manageable amount of time. 

The gas distributions in the p-T plane at this end state are illustrated in Fig. 5. At densities less 
than ~ 5 X 10~^^g cm~^, the distributions for the 32^ (top left) and 64^ (bottom left) simulations 
are identical in form. At ~ 5 x 10~^^ g cm ~^ , it is the lower mass-resolution simulation that 
produces the greatest amount of cool, dense gas. 

Evolving these sets of data while enforcing a minimum temperature sufficient to pressure 
support structures with NgpH or fewer particles does little to alter these distributions (Fig. 5, 
panels to the right) outside those regions below T^m and below the line with a slope of Vs which 
extends from the point of initial density and temperature of the gas. This line corresponds to the 
path a uniform gas would take due to the imposition of a minimum temperature. 

Recall from Sec. 3.1 that individual structures in high resolution simulations have, at early 
times, less shocked material in their cores than is produced at lower mass resolution. This difference 
is eliminated after the initial merging period on account of the introduction of new structures falling 
into the lower-resolution cluster. These structures have material with low degrees of shocking. 
Compare the lower panels of Fig. 1 with Fig. 6 (left panels) which plots the p-T distribution 
for a cluster from a 32^ and 64^ simulation at late times (note that the figures have different 
scales). At early times, there is a population of less shocked gas in the high resolution runs 
as described earlier. This leads to the spread in temperature at high densities. Although the 
distribution of substructures differs among the simulations of differing resolution, the presence 
of substructures at the termination of the simulations leads to a similar spread in temperatures, 
regardless of resolution. Taken as singular objects, the bulk properties of individual clusters do 
not vary substantially between simulations of differing resolutions. The degree of shocking, as 
parameterized by the entropy factor < Tp~^ >m (a mass- weighted mean) is found to be essentially 
constant for the cluster simulated at either resolution. It is (1.12 it 0.02) x 10^^ K cm^g"^/^ for 
both the 32'^ and 64^ simulations. With the imposition of a minimum temperature, the entropy is 
essentially unchanged at (1.13 ± 0.03) x lO^^K cm^g^^/^ The trend for the mean temperature of 
structures to increase with increasing resolution is maintained, with the gas temperature rising from 
(2.9 ± 0.1) X 10^ K to (3.2 ± 0.1) x 10^ K between the 32^ and 64^ simulations. This is understood 
to be a by-product of the deeper potential wells in the cores of structures as described earlier. 

The density distribution confirms the similarity of the states of the gas in the simulations 
of differing resolution (Fig. 3, solid line). In the absence of a minimum temperature (Fig. 3, 
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left panels), the gas density distribution varies marginally with resolution with a decrease in the 
amount of gas at a density of ~ 10~^^ g cm ~^ compensated by an increase at the higher density 
of ^ 10^^^ g cm as the spatial resolution is doubled. The addition of a minimum temperature 
skews the distribution even more significantly to lower densities for the lower-resolution simulation. 
The shift in the simulations without a minimum temperature can be explained as a result of the 
greater abundance of substructure in the dark matter distribution which locally compresses and 
heats the gas. The gas at the highest densities in the cores of structures {p > 10~^^g cm~^), is 
not significantly affected. 

Finally, a useful diagnostic is the mean density profile (Fig. 2) of the clusters at the final epoch, 
a = ttf. It confirms that the presence of a minimum temperature does not change the morphology 
of the clusters. The mean profiles are calculated from a sample of clusters selected by their size, 
to ensure that the profiles are resolved. The inner cores are excluded from the calculation of the 
means, since they are not resolved. 

The simulations of OV97 display a marked tendency for the gas to become enriched in the 
cores of clusters, illustrated in Fig. 9 of OV97. This is contrary to the results found by others 
(Evrard 1990; Thomas & Couchman 1992; Cen & Ostriker 1993; Kang et al. 1994; Metzler & 
Evrard 1994; Pearce et al. 1994; Navarro et al. 1995; Anninos & Norman 1996; Lubin et al. 1996; 
Pildis et al. 1996). Indeed, it is most commonly found that the gas becomes more dispersed than 
the dark matter. This result is supported by the simulations we performed. The baryon fraction 
parameter, T, which is the ratio of the baryon density to total mass density normalized by the 
global mean ratio, is plotted in Fig. 7 for the clusters. The profiles given axe the mean profiles for 
the same sample of clusters for which the mean density profiles have been calculated. It is also clear 
from this figure that the presence of a minimum temperature has no effect on the baryon fraction 
distribution. It also shows clearly that resolution does have an effect in the cores. 
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4. DISCUSSION 

It has been suggested (Owen &; Villumsen 1997) that convergence with increasing mass res- 
olution in numerical simulations of hierarchical structure formation requires the imposition of a 
minimum temperature in order that gas is pressure supported in structures with masses similar to 
the effective resolution of the simulation. Our results do not support this contention, either at early 
epochs after only a few mergers, or at later epochs after a significant number of mergers have oc- 
curred and the minimum temperature has dropped well below the virial temperature of the typical 
structures. Our results indicate that the problem of the gas becoming over-shocked in simulations 
in which the Jeans Mass is not resolved is not as severe as suggested by OV97. The larger degree 
of shocking of the bulk of the gas noted by OV97 as resolution is lowered is not reproduced at any 
epoch of the simulations presented here. Consequently, the imposition of a minimum temperature 
is not necessary to achieve convergence. 

We do observe a phenomenon which is similar to the decrease in shocking with increasing 
resolution as described in OV97. At an early epoch, isolated halos are shown to possess a large 
amount of substructure in their high resolution simulations that is not present in their low resolution 
counterparts. This substructure contains cool, dense, weakly shocked gas. However, the presence 
of this population of cool and dense gas in the high resolution simulations is offset by an increase 
in the temperature of the halo gas surrounding these substructures. The net effect is not to 
lower the mean temperature of the halo but to marginally increase it (10% for each doubling in 
resolution) as well as to increase the entropy of the halo by an equivalent fraction. At a later epoch, 
after subsequent accumulation of material, the halos in the lower resolution simulations develop a 
comparable population of cooler, unshocked gas within substructure. 

Further attempts were made to reproduce the results of OV97 by modifying our code. Since it 
was believed that the gravity calculations for the simulations of OV97 were done using a particle- 
mesh method, we disabled the particle-particle contributions of the gravitational force as well as 
the refinements in the HYDRA code. This was unsuccessful in reproducing their results. 

The results of imposing a minimum temperature sufficient to pressure-support the resolution- 
limited gas structures were examined for adverse effects. At the final epoch, the gas density profiles 
of the clusters were not altered. The degree of shocking increased by less than 5%. This increase 
may be attributed primarily to an increase in temperature and decrease in density of the gas lying 
near the cores of structures. The bulk of the gas, which lies in the halo, was unaffected. 

The concentration of baryonic material in the cores of structures reported by OV97 is also 
not observed. Instead, an anti-bias is found, in accordance with the substantial majority of nu- 
merical simulations reported in the literature as noted previously. Pearce et al. (1994) explains 
the phenomenon as a result of the merging process of gas-dark-matter halo systems in which gas 
is shocked, permanently removing the energy from the dark matter component. For the gas to be 
more biased in clusters, without an energy-loss mechanism operating such as cooling, the gas must 
transfer energy to the dark matter component. 
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It is unclear what causes the differences between our results and those of OV97. The most 
striking remaining difference between our investigations and those of OV97, which we have not ex- 
plored, is the difference in dimensionality of the numerical simulations. OV97 used two-dimensional 
simulations, ours were three-dimensional. Structures in two dimensions are rods, which have es- 
sentially the same cross sections per mass as spheres, but one less degree of freedom in which to 
avoid collisions. This may strongly affect the efficiency of shock heating, with more head-on col- 
lisions per merging event leading to more violently shocked gas. However, naive interpretation of 
this observation would suggest that the gas halos should be more extended, and more anti-biased 
relative to the dark matter, contrary to their results. 

We conclude that resolution of the Jeans mass of the smallest structures is consequently not 
important in studies of the gas in structures after merging has occurred. It may be important for 
the first objects formed, but for any study other than the statistics of position and number, these 
first objects should be ignored since they are, by definition, poorly resolved. Clearly, this result 
may not hold for simulations in which there is cooling and/or star formation feedback, but these 
will have an effect only in the densest parts of the clusters, and should not affect global properties. 
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Table 1. Simulation properties 



# of particles 


e 


Resolution gas dark {h 


-^kpc) 


Evolved to a = 0.17a/ 


32^ 32^ 32^ 


80 


643 543 


40 


1283 1283 1283 


20 


Evolved to a = 0.22a/ and 


«/ 


323 323 


80 


643 643 643 


40 



Table 2. The measure of shocking of an individual cluster 





<T/p"J' >M (K cm2g-V3) 


Resolution 


Unsampled Sampled and recalculated pg 


323 


5.5 X 10^4 5.5 X 10^4 


643 


5.9 X 10^4 1.6 X 10^4 


1283 


6.5 X 10^4 0.5 X 10^4 
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Fig. 1. — Temperature vs. gas density at a timescale factor, a = 0.17aj. The upper set of panels are 
for all the gas while the lower set is for a selected cluster. From left to right, the panels correspond 
to simulations with increasing resolutions of 32^, 64"^, and 128^ baryon particles. The simulation 
data were resampled before analysis to the same mass resolution as the 32^ simulation. Contour 
lines in the upper panels are of constant mass surface density in p — T space. In the lower set of 
panels, the data points are for the sample of particles in the selected cluster. 
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32^ 64^ 




Fig. 2. — The mean dark matter and gas density profiles. The panel to the left corresponds to 
the low resolution set of runs, while the panel on the right is for the high resolution set. Profiles 
are given for clusters formed in the absence of a minimum temperature and for clusters formed in 
a simulation with a minimum temperature scaled appropriately for the mean gas density at each 
epoch. 
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Fig. 3. — The distribution of gas density at early and late epochs. The upper set of panels is for 
the 32^ set of data while the lower set is for the 64^ set. From left to right, the panels correspond 
to no minimum temperature and a minimum temperature scaled appropriately for the epoch. 
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Fig. 4. — Temperature vs. gas density for all the gas at an expansion factor of a = 0.22a/. The 

circles in the right-hand panels indicate the starting density and temperature. For the left-hand 
panels, these points are off the figure at T = 100 K . The lines projecting to the lower left of each 
circle correspond to the path of taken by a uniform gas expanding and cooling adiabatically (lower 
line) and to the path of gas maintaining Tmin (upper line). The panels are arranged as in Fig. 3. 
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Fig. 6. — Temperature vs. gas density for the gas of a single cluster at an expansion factor, a = af. 
The panels are as described in Fig. 4. 
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Fig. 7. — The radial profile of the mean baryon fraction parameter, T, for the clusters formed at 
different resolutions in the presence or absence of a minimum temperature. 



